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The Hawking effect can be rederived in terms of two-point functions and in such a way that 
it makes it possible to estimate, within the conventional semiclassical theory, the contribution of 
ultrashort distances at Z"*" to the Planckian spectrum. The analysis shows that, for Schwarzschild 
astrophysical black holes, the Hawking radiation (for both bosons and fermions) is very robust up 
to very high frequencies (typically two orders above Hawking's temperature). Below this scale, the 
' contribution of ultrashort distances to the spectrum is negligible. We argue, using a simple model 

with modified two-point functions, that the above result seems to have a general validity and that it 
is related to the observer independence of the short-distance behavior of the corresponding two-point 
^ ' function. The above suggests that only at high emission frequencies could an underlying quantum 

Q , theory of gravity potentially predict significant deviations from Hawking's semiclassical result. 
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■ I. INTRODUCTION 

in 
m 



Semiclassical gravity predicts the radiation of quanta by black holes [Tl, Q ■ The emission rate is given by the product 
of the Planckian factor times the grey-body coefficient Timp{w) 



O. dNi„^p{w) _ 1 1 

> ] 

JL . where k, Qh and are the surface gravity, angular velocity and the electric potential of the black hole horizon. 
The signs ± in the denominator account for the Bose or Fermi statistics, and m, p and q are the corresponding axial 
, angular momentum, helicity, and charge of the radiated particle. 
;V ■ The deep connection of this result with thermodynamics [Sj] and, in particular, with a generalized second law [3| 
• 5^ I strongly supports its robustness [1, H, [3|- However, as stressed in a crucial ingredient in deriving Hawking 
' radiation via semiclassical gravity is the fact that any emitted quanta, even those with very low frequency at future 
^ infinity, will suffer a divergent blueshift when propagated backwards in time and measured by a freely falling observer. 
Also, in the derivation of Fredenhagen and Haag @, the role of the short-distance behavior of the two-point function 
is fundamental. All derivations seem to invoke Planck-scale physics. The exponential blueshift effect of the horizon 
of the black hole could thus be regarded as a magnifying glass that makes visible the ultrashort-distance physics 
to external observers. According to this reasoning the microscopic structure offered by string theory (or any other 
underlying theory) could leave some imprint or signal in the emission rate. However, the results of string theory 
seem to agree with Hawking's prediction. For the emission of low-energy quanta (with wavelength large compared 
to the black hole radius), and for some particular near-extremal charged black holes, the prediction of string theory 
[13, [iH coincides with the rate ([T]). This agreement is complete, despite the fact that the two calculations are very 
different. For instance, whereas in semiclassical gravity one can naturally split the emission rate into two factors 
(pure Planckian black-body term and grey-body factor), in the D-brane derivation one gets directly the final answer 
without the above mentioned splitting. 

While the calculation of string theory requires low frequencies for the emitted particles, in the arena of semiclassical 
gravity the result is valid for all wavelengths, even those smaller than the size of the black hole. The thermodynamic 
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picture strongly suggests the robustness of Hawking's prediction and its interpretation as a low-energy effect, not 
affected by the particular underlying theory of quantum gravity (see also [l^l), and expected to be valid for a large 
range of frequencies. However, from the perspective of quantum field theory in curved spacetime, it is unclear how to 
introduce a cutoff in the scheme (parameterizing our ignorance on trans-Planckian physics) in such a way that, for low- 
energy emitted quanta, the decay rate ([T]) is kept unaltered. The aim of this work is to study this issue in some detail^. 

In section |lT] we will review the standard derivation of black hole radiation emphasizing the role of ultrahigh 
frequencies to get the Planckian spectrum. In section IIIII we rederive the Hawking effect in terms of two-point 
functions, instead of Bogolubov transformations (for a general reference see [I^l), for both massless scalar and spin- 
1/2 fields. The new derivation of the black hole decay rate offers an explicit way to evaluate the contribution of 
ultrashort (Planck-scale) distances to the thermal Hawking spectrum. This is the subject of section IV. In section 
V we present a simple model, where the two-point functions are deformed with a Planck-length parameter, to show 
how the previous results emerge in this new scenario and support their robustness. We point out that a generalized 
Hadamard condition plays a fundamental role to keep unaltered the bulk of the Hawking effect. Finally, in section 
VI, we summarize our conclusions and make some speculative comments. In the appendices we give details of some 
calculations used in the body of the text. 

II. BOGOLUBOV COEFFICIENTS AND BLACK HOLE RADIANCE 

Let us consider the formation process of a Schwarzschild black hole, as depicted in Fig.l, and a massless real scalar 
field (j) propagating in this background. The equation of motion obeyed by the field is Ocf) = and the Klein-Gordon 
scalar product is given by 

(01,02) - / rfS^(0i5^0; - 0;9^0i) , (2) 

where S is a suitable "initial data" hypersurface. A natural choice for S is the past null infinity /~ and therefore one 
can express the field in a set of modes u™(a:), which have positive frequency in 

1 \ ^ / in in , in] in^\ /o\ 

9^ 2^{a^ u, 'u, ) . (3) 

i 

Alternatively, we can choose E as E = /+ U , where is the future null infinity and is the event horizon. 



r=0 




Figure 1: Penrose diagram of a collapsing body producing a Schwarzschild black hole. 



We note that this sort of problem has aheady been addressed in the context of acoustic black holes by modifying the dispersion relation 
of the wave equation obeyed by sonic disturbances |l4l llSlI . This is naturally justified as an effect of the atomic microscopic structure 
of the fluid, and requires a breakdown of Lorentz invariance. The rest frame of the atoms of the fluid plays a privileged role. In this 
paper we follow an alternative route. 
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According to this we can then expand the field in an orthonormal set of modes u°^*{x), which have positive frequency 
with respect to the inertial time at /+ and have zero Cauchy data in together with a set of modes with 
null outgoing component at /+. Therefore we can write 

^ = ^(ar*<"* + ar*^ur**) + (ar^^r* + ar*^«r**) ■ (4) 

i 

The particular choice of modes m*"* does not affect the computation of particle production at /+ , so we leave them 
unspecified. 

The modes M°"*(a;) can be expressed in terms of the basis u™ 

uf^\x) = «.«"r(^) + > (5) 

i 

where the coefficients aji and pji are the so-called Bogolubov coefficients and are given by the scalar products 

a., = {uT\uf) , - -«"*, uf*) . (6) 
The above expansion leads to an analogous relation for the creation and annihilation operators: 

„out \ '' f„* „in Q* „i"t\ l'7\ 

When the coefficients /3y do not vanish the vacuum states \in) and \out), defined as a*"|m) — and a°'^*\out) = 0, 
do not coincide and, as a consequence, the number of particles measured in the i^^ mode by an "out" observer, 
7V°"* — /i~^a°"*^a°"*, in the state \in) is given by 

H7vr*H = ^iAfcr (8) 
fc 

Let us now briefly summarize the main steps of Hawking's derivation. Assuming for simplicity that the background 
is spherically symmetric we can choose the following basis for the ingoing and outgoing modes 

v^:iji — = — Yr{9,^). (9) 

1 p—iwu 

uirji^ ^ -= — Yr{e.4>) . (10) 

V47rw r 

Here Y^(6,4>) are the spherical harmonics. One can evaluate the coefficients Pwim,w'i'm' according to previous 
expressions by making the convenient choice T, = I~ 

f3wim,w'i'm' = i J dvr'^dQ{u'^l^dvU^^q,„^, — u™,i,j^,dvu'^l^) . (11) 

The angular integration is straightforward and leads to delta functions Sw S„i^^m' for the P coefficients. The relevant 
point is to realize that the coefficients can be evaluated and have a unique answer, which turns out to be independent 
of the details of the collapse, if represents a late-time wave-packet mode (i.e., centered around an instant u with 
u +00 along /+). When these modes are propagated backwards in time they are largely blueshifted when they 
approach the event horizon. After passing through the collapsing body they are scattered to /~ in a very small 
interval just before vh- To know how they behave on (as needed to evaluate the scalar product with m^;^) one 
can apply the geometrical optics approximation since the effective frequency, as measured by freely falling observers, 
is very large. The (late-time) mode which is of the form pH)) at /+, evolves and arrives at /~ with the form 

uifji ^ yne^cPMvH - v) , (12) 

where t; (w) is the transmission coefficient for the Schwarzschild metric and the relation between null inertial coordi- 
nates u at and v at /~ is typically given by the logarithmic term 

u = Vh — \n k\vh — v\ , (13) 
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where, for the Schwarzschild black hole, k = 1/4M and vh represents the location of the null ray that will form the 
event horizon iJ+ (see Fig. 1). One has then all ingredients to work out the (late-time) Bogolubov coefficients 

B I n, , — Hilll— M^^^W^ I rf,,p-™(fH-K"MnK|tiH-t>|)-itu'-u r c 

277 y W J_^ 



They can be evaluated explicitly 

Pwlm,w'l'm' = -^—L-HnL.I— " , i,„T (l + K~'wz)(5n"^m-m' , (15) 

ZTTK V W (— K '-W'l + 6)'- + '^ 

where we have introduced a negative real part (— e) into the exponent of (|14p to ensure convergence of the corresponding 
integrals. To compute the particle production at one has to evaluate the integral (from now on in this section we 
shall omit, for simplicity, the subscripts I, m) 

+00 

dw'Piuiw'P^^^' ■ (16) 

Jo 

The integration in w' reduces to 

dw' 







from which we finally get 



2TTKe-""' '"'5{wi-W2) , (17) 



dw'p^,^,pi^^, = J!;'}^^l!}'_^ 5{w, - W2) , (18) 



where the coefficient in front of S{wi — W2) represents a steady thermal flow of radiation of frequency w = wi 

^^54^ - ^{^nlKTM = J- , ^'1") , (19) 
and the grey-body factor is given by Ti{w) = |t/(w)p. For a generic collapse the result leads to formula H]). 

It is important to remark at this point that a basic step to exactly obtain the Planckian spectrum is (|17|) . which 
crucially requires an unbounded integration in all frequencies w' . In fact, if we introduce an ultraviolet cutoff A for 
w' we should replace (fT7|) by 

dw^^-^-^n,^^ln(-^-^w'-^e)^K-^W2ilni^-'w'-^e) ^ g"''''" '""^ 27r,5^ [k"! (lUi - W2)] , (20) 

w 

where we have defined 

sin r^SlVi^" 

Sa[K-\wi ~ W2)] - ^-jj-^ f (21) 

7rK~^(wi — W2) 

ln[K-iA] ^ ' 

Note that in the limit as a goes to zero 8^ turns into Dirac's delta function and we recover (|17|) . The new expression 
is, however, qualitatively different from the previous one. To evaluate the new emission rate requires making use of 
normalized wave-packet modes. Introducing the standard ones [l| 







1 /-U + iJc 

— / p2mwn e out 

— I aw e ^wim I 

V e Jje 



'^'jnim — ~~i= I dw e ^ , (23) 



where j > and n are integers, representing wave-packets peaked around the retarded time u„ = 27rn/e and centered, 
with width e, around the frequency Wj = (j -|- l/2)e; and, accordingly, defining 

/3,„,^- = ^ / dw e^™"/' p^^, , (24) 
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the emission rate results (see appendix A): 



1 [(2^_^^)^ 



From this expression^ we see that the rate of emitted particles depends on the retarded time w„ = 27m/e and decays 
with time for any small but nonzero value of cr = l/ln[A/K]. Only when A goes to infinity (no high frequency cutoff) 
do we recover the steady thermal flux of radiation. 

In conclusion, the above discussion shows that the radiation is now time-dependent and decays for all finite values 
of A. The decay in time would also occur at low frequencies, where string theory agrees with Hawking's prediction. 
Therefore, as expected from conventional arguments, the mathematical role of the ultrahigh frequencies is very 
important for the late-time behavior. Nevertheless, since they only enter as virtual quanta, their actual status is 
unclear [l6j . A derivation of the Hawking effect based on quantities defined on the asymptotically flat region, where 
physical observations are made, would be preferable. This turns out to be possible if, instead of working with 
Bogolubov coefficients, one uses two-point functions. They are defined in the /+ region where a Planck-length cutoff 
in "distances" can be naturally introduced. This is the task of next sections. 



III. TWO-POINT FUNCTIONS AND BLACK HOLE RADIANCE 



This section will be devoted to rederive Hawking radiation by means of two-point functions. Intuitively the idea is 
simple. In the conventional analysis in terms of Bogolubov coefficients, we first perform the integration in distances 
(to compute the scalar product required for the f) coefficients) and leave to the end the integration in frequencies 
w' . In contrast, we can invert the order and perform first the integration in frequencies (which naturally leads to 
introduce the two-point function of the matter field) and perform the integration in distances at the end. 

Let us rewrite the basic expression ([5]), or more precisely, the expectation values of the operator iV^""* = 

tr^ a°"*|^a°"*, as follows 

k k 

= E (//^i<"*(^i)^M«r(^i)) (£ds>°f*(x2)a.ur(^2)) . (26) 

If we now consider the sum in modes before making the integrals of the two scalar products, and take into account 
that 

(zn|</.(xi)0(x2)|m) = hJ2ut{xi)uT*{x2) , (27) 

k 

we obtain a simple expression for the particle production number in terms of the two-point function 

{in\N,{t\^n) = [ dSi'dS^Kr*(2^i)^,JKr(x2)5,](zn|0(a;i)0(a;2)|*n) . (28) 

In the above expression the two-point function should be then interpreted in the distributional sense. The "ie- 
prescription" (see ea. psp below) is therefore assumed for the two-point distribution {in\(j){xi)4>{x2)\in) and it verifies 
the Hadamard condition"^ i^, Alternatively, taking into account the trivial identity {out\a°^*l_^a°^*^\out) = we 
can rewrite the above expression as [l8| 

{in\N°X\in) = h-' f dI]5'dS^[<'*(xi)a^]Kf *(a;2)a,](*n| : c^{x,)<l>{x2) : \in) , (29) 

where {in\ : (j){xi)(f>{x2) '■ \in) = {in\(f>{xi)(l){x2)\in) — {out\(j)(xi)(j){x2)\out) . Now the Hadamard condition for both 
\in) and \out) states ensures that {in\ : (p{xi)(j)(x2) ■ \in) is a smooth function. 



^ We note that the oscillatory behavior in 1251 1 is an artifact of the particular way we have introduced the cutoff. If the cutoff is introduced 
in a different way, see appendix A, the oscillatory term disappears but the decay with time is maintained as ^ e~^^'^'^"^'^~"H)('^i^'^ Z'^)] , 

^ The two-point distribution should have, for all physical states, a short-distance structure similar to that of the ordinary vacuum state 
in Minkowski space: {2it)~^{(j + 2iet + e'^)~^ , where a{x\, X2) is the squared geodesic distance. 
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A. Thermal spectrum for a scalar field 

Let us now apply this scheme in the formation process of a spherically symmetric black hole and restrict the "out" 
region to J+. The "in" region is, as usual, defined by /~. At /+ we can consider the conventional radial plane- wave 
modes 

uir^{t,r,e,cl>)\i^ ~ uZ\u)^^^^^ , (30) 

where = ^/^=- We shall now evaluate the matrix coefficients (m|-/V°"*|m) where zi,2 = (wi,2, ?i,2, "11,2)- 

Taking as the initial value hypersurface /~ and integrating by parts we obtain 



out^,oul^ ^ 

W2 ^ 



ri r2 

The two-point function above can be now expanded at /~ as 
{in\(t>{ 



dy^dy2{in\(j){xi)(j){x2)\in) . (31) 



xi)4)(x2)\m) = h dw} . — ■j^=— -. (32) 

l,m 

Recall that the radial part of the late-time "out" modes w°"m' when they are propagated backward in time and reach 
/~ , takes the form 



e 



-iwu{v) 



uT\i- ~ ti{w) _ e{vH - v) (33) 
where u{v) ~ vh — In. k{vh — v). Performing now angular integrations and taking into account that 

-iw(vi—V2) ^ 2 



(9 (9 



/ dw—^ - - — -2 , (34) 

^0 47r'u; 47r [vi — V2 — te) 



we get 



tl (W-\)tr iWt) I""" p-lWlUlVl)+tW2U(_V2) 

{in\N?X\in) = - \\^Jr_J' / dv,dv2—, r^Si,iJm,m2 , (35) 

where the limit e ^ 0+ is understood. Alternatively, since we are interested in quantities measured at S = /+, we 
could use this latter hypersurface to carry out the calculations. In this case, the expression for the particle production 
rate becomes 

(^nKtl-) = r d..c^.2 ,,^^^^yt^^, e----"-^---"-.,,^^.^. , (36) 

A-K^y/wrm J -00 [v(ui) - v{u2) - leY 

{in\NtX\in) = 'iii-'^ / du,du2j^^, -^KiM,m2 ■ (37) 

47r^^wiW2 J-00 [smh f - 7i2 - 2e)J2 

This last expression is more convenient for computational purposes. Since the function in the integral depends only 
on the difference z = U2 — u\, the integral in U2 + ui can be performed immediately and leads to a delta function in 
frequencies 



and leads to 



. = _ ^^.(-0y^^^.-^2) p,.e--4-^ , (38) 

2TTy/wrw^ i-cx> [smhf (2-ie)]2 



tlAwi)tl{n,2)5{w^ - W2) /■+°° {%?Kl25m^r 

-( 

2 > 

Performing the integration in z we recover the Planckian spectrum and the particle production rate 



(zn|iV:-|zn) = J^M^. (39) 
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This derivation of black hole radiation is somewhat parallel to the one given in Q . The emphasis is in the two-point 
function of the quantum state, instead of the usual treatment in terms of Bogolubov transformations. It is worth 
noting that ([35| displays an apparent sensitivity to ultrashort distances due to the highly oscillatory behavior of the 
modes in a small region before vr- A similar conclusion can be obtained from (|38p when z — > 0. The sensitivity to 
short distances is, however, less apparent if we repeat the above calculations using the expression ([29]) instead of (|28)) . 
In this case, we find 



ti,{wi)ti{w2) r" , . 

- . o , — =^ / dvidv2e 
1 TijMltM 



{Vi-V2)^ [u{vi) - U{V2)Y 



-iwiu{vi)-\-iw2u{v2) 



(40) 



where we have dropped the ie terms since they are now redundant. Note that the short-distance divergence of 



— 1)2)^ in pO]) is exactly cancelled by 



7 ("2) 



-/ — , — 5 — for any smooth choice of the function u(v). This cancellation is 
a consequence of the Hadamard condition that verify both "in" and "out" vacuum states. It is also important to remark 
that the above formula exhibits the absence of particle production under conformal-type (Mobius) transformations 



(41) 



where ah — cd — 1.^ 

If the calculation is performed using E = /+ , one finds 



th{wi)t*^{w2) 

4:W^y/WlW2 

fiK)t("2) 
{v{ui) - viu2)y 



[Ui - U2 



(42) 



which leads to 



(zn|iV°tKn) 



tliiwi)ti {W2)S(WI -W2) 



2tT^WiW2 



-t-00 



dze 



(sinh ^z)2 

The integral in distances z also leads, as expected, to the Hawking formula^ 



''l+"'2 . 



(43) 



(inlN, 



wlm I 



\tiH\' 



+00 



i-2? 



(sinh|z)2 



(44) 



B. Thermal spectrum for a s — 1/2 field 

In this subsection we shall extend the analysis of the scalar field to a fermionic s = 1/2 field. For simplicity we take 
a massless Dirac field, obeying the wave equation 

7^V^7/' - , (45) 



* This leads, immediately, to the expected result that there is no particle production under Lorentz transformations. 
For Kerr-Newman black holes the calculation is similar up to a shift in the wave function e~""^, which should be now replaced by 
g— 2(jD-mnij— (j*ff )z^ g^j^ effect of wave propagation through the corresponding potential barrier. 



8 



where 7^ = V^j"' are the curved space counterparts of the Dirac matrices 7° (see appendix B for calculations omitted 
in this section). The Klein-Gordon scalar product ^ is now replaced by 



= / dS'^^i7^V2 . (46) 
Therefore the expression for the expectation values (pS)) is replaced by 

{in\N^X\'n) = fi-i / dS5'dS^[<"*(a:2)7.]6[7M<"*(^i)]"(*«l'/5a(a:i)V'''(2:2)Nn) • (47) 
is 

At J"*" we can consider the normalized radial plane- wave modes^ 

-iwu f \ 

<^.^.^^^^^M..-^{^;^ , (48) 

where 7?(r)™^ are two-component spinor harmonics (see appendix B). Note that the angular momentum quantum 
number j is uniquely determined by the relation Kj = ±(j -f 1/2). The above modes, when propagated backwards in 
time and reach , turn into 



out 



where the last term y^du{v) / dv appears due to the fermionic character of the field^. Proceeding as in the bosonic 
case we can expand the two-point function as 

{in\4>a{xi)'ip''{x2)\in) = h^v]^';^{xi)vj!'''' {X2) (50) 

A; 

where ti™ are negative-energy solutions which in take the form 

< - <%„.(^,^,M)I.- - ^ ( ) , (51) 

Performing first the angular integrations and taking into account the orthonormality relations of the spinor harmonics 
vi^)^/ 1 the above formulas get simplified and become 



X / dvidv2\l— j-r- (52) 

00 



dv dv {vi ~ V2 ~ it) 
As in the bosonic case, we rewrite this expression as an integral over 

/■ l^routl- N .iK,,(Wl)t* (W2) ^ 



X 



00 



— 00 



dMidM2e-™i"i+™^"^— 7-^-^^^ -T (53) 

smh[^(ui — U2 — ie)J 



We can also split the integral in a product of a function dependent on U2 + u\ and another function which depends 
on z = M2 — Ml. The former leads to a delta function in frequencies and we are left with 

= -^^«:,i('«iK,Jw2)(5™^^™^/«^.^«^.^(5(«;i - W2) X 

+ 00 IK\ 

dze-^^^ ^ ; 

smh[f(z-ze)] 



® On physical grounds we should use left-handed spinors u°j^^-^_ = ^ ('''JJ^"*^ . | „ . ~ ""w^-Ik |m '^'^^ final result is not changed. See 
appendix B. 

The vierbein fields needed to properly write the field equation in a curved spacetime have been trivially fixed in the asymptotic flat 
regions, so its transformation law under change of coordinates is then translated to the spinor itself (see also appendix B). 
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Performing now the integration in z 



—I 
2^ 



+ 00 



dze 



(f) 



sinh[f (z - ie)] e^™" ' + 1 



(55) 



we recover the Planckian spectrum, with the Dirac-Fermi statistics, and the corresponding particle production rate 



Analogously as in the bosonic case, if we use the normal-ordering prescription instead of the ie one we find 



47r2 

-iwiu{vi)-\'iw2u{v2} 



,. |,.out,. X tK,,{wi)t^ {W2)* r» I du{vi) du{v2) 

{i-nWiti^Vn) = -I— -^r^ / dvidv2\l — — T-^ X 



dv dv 

du(vi) du{v2) 
dv dv 



[V1-V2] [u{vi) - U{V2)] 



(56) 



(57) 



Jl '^J2 Jl J2 



Changing the mtegration surface to /"^ we get 



^l(Wl)tK,, {W2)* 



47r2 



duidu2 



1+ 



-iwiU-i+iW2U2 



dvi^ui) dv{u2) 
du du 



dv{ui) dv{u2) 
du du 

1 



(58) 



[v{ui) ~ V{U2)] [U1-U2] 



^'^Jl '^32 ^^31 ^32 



Note again that the short-distance divergence of the two-point function of the "out" state ^^^^^^j is exactly cancelled 

by the corresponding one of the "in" state [y[ui)-v{u2)] ' ^'^^^^ both vacua are Hadamard states. After some algebra 
we get 



\t (w)\^ r^°° 



(f) 



1 



and taking into account 



+00 

— / dze"™^ 
27r 



(f) 



sinh(-|z) z 



sinh(^z) z 



+ 1 



(59) 



(60) 



we newly recover the fermionic thermal spectrum. 



IV. SHORT-DISTANCE CONTRIBUTION TO THE PLANCKIAN SPECTRUM 

A. Bosons 

We have seen in the previous section that it is possible to rederive the Hawking effect in terms of two-point functions. 
Either via expressions (|^ . ([55)1 or, equivalently, via expressions (P5t .p ^ . Both prescriptions are equivalent and lead 
to the Planckian spectrum modulated by grey-body factors. The advantage of the final expression ([H]) is that it offers 
an explicit evaluation of the contribution of distances to the Planckian spectrum. To be more explicit, the integral* 



(wK ^ , an) — — 



2lTW 



dze 



(sinh |z)2 



(61) 



We have intentionally omitted the grey-body factors in 1611 1 because they are irrelevant for the discussion of this section. 
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can be interpreted as the contribution of short-distances z € [—a, a] to the (bosonic) thermal spectrum when a is 
close to the Planck length Ip. One could, alternatively, be tempted to propose, according to the integral 

' 2ttwJ_^ [sinh|(z-ie)]2 ^ ' 

as a legitimate expression to account for the short-distance contributions. However, this interpretation is not physically 
sound. In the absence of a black hole, when there is no radiation at all, the above expression becomes 

-1 /■+" 1 

dze"™"- . (63) 



27rw J_a — 

For a — > +00 this expression vanishes, as expected due to the absence of radiation. However, for finite a it is 
non- vanishing. In contrast, the proposed expression (j61[) . does not suffer from this weird behavior, due to the 
presence of the second term. 



In conclusion, the calculation of black hole radiation using the prescription ([29)1 offers the possibility to re-evaluate 
Hawking radiation by removing the range of distances where physics can be dominated by an underlying theory 
beyond field theory. We shall now work out explicitly the short-distance contribution to Hawking radiation to see 
whether it is fundamental or not in order to obtain the thermal spectrum. The integral (j6ip can be worked out 
analytically 

/^(u.«-\ aK) = -^l!^ _ J:L{e"-(F[l, ^iwK-\ 1 - iwK~\e""'] 



Attw 



^F[l,-iwK ^,1-iwK ^ , e'^'^])} + — cos(q;u;) 

ZTTaW 



aK-, — 2 

(gQK _ 2) 



(64) 



where F is a hypergeomctric function and Si{x) — dt^^. To get some insight about the properties of this formula, 
we find useful to expand it in powers of wk~^ and an. The expansion in an assumes that the microscopic length 
scale a Ip is much smaller than the typical emission wavelength ~ of the black hole, whose (macroscopic) 
temperature is Th — Kjl-K. For a Solar-mass black hole aK ~ 10~^° and for a primordial black hole of 10^'^ g 
aK ~ 10^2^. The expansion in wkT^ means that we are looking at frequencies below the typical emission frequency, 
wtypicai ~ Th, of the black hole. The result is as follows 



I^(wK^^,aK) = ( -^an —(aK)^ + 0\(aK)'^]\ — 

^ ' ' \127T 720n ' ' w 



- i- 



)-{aKf + 0\{aKf\\ - -f {p\{aKf\) {-f + . . . (65) 
27r J K K 



\72tt 

From this expansion we conclude that the contribution of short distances to the spectrum is completely negligible in 
the very low energy regime w/k <ti 1 since 



linr //-"^ '""^ = — «1. (66) 



Moreover, due to the smallness of an, we find that {wtypicaii^ ^,Ka) can be well approximated by (|65p even for 
frequencies close to the typical emission frequency, which leads to 

I^iw,yp..a^K-\aK) _,03^^^<^ (g^) 



Again, since an <^ 1, we find a negligible contribution at Wtypicai ^ Th- To be precise, for a Schwarzschild black 
hole of three solar masses, when a is around the Planck length = 1.6 x 10^'^^cm, the relative contribution to the 

Planckian distribution '"^^l^ is, for w ~ wtypicai, of order 10"'^^%. For primordial black holes, M ^ 10^^ g, 

the relative contribution is still insignificant: 10^^^%. Even more, using the expansion (j65p we easily get 

I^{wK-\aK) _ c,K(e27r^,K-i _ 
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and we find that, for a black hole of three solar masses, we need to look at the high frequency region, w/wtypicai ~ 96, 
to find that the contribution of Planck distances I^{wk~^, Ipn) is of order of the total spectrum itself® . For primordial 
black holes we find w / Wtypicai ~ 52. The same numerical estimates can be found using the exact analytical expressions. 

We can also naturally ask about the contribution to the spectrum of large distances. This question is immediately 
answered using our analytical expression (|64|) . The contribution of distances up to a = 20rg, where Vg is the 
gravitational radius, represents 90% of the thermal peak at Wtypicai- For a — 200rg we obtain 99.7% and for 
a = 2 X lO'^rg the percentage is around 99.99998% . 

Summarizing, we have provided a quantitative estimate of how much of Hawking radiation is actually due to 
Planckian distances. It turns out that the contribution of ultrashort distances is negligible and thermal radiation 
is very robust up to frequencies of order 96Th (for Schwarzschild black holes of three solar masses) or 52Th (for 
primordial black holes). In parallel and dual to this, the contribution of large distances is also insignificant. 

It is interesting to repeat the same calculations with the ie-prescription. As we have already stressed with this 
prescription one cannot expect a meaningful result. The outcome is completely different. The contribution of distances 
in the interval z £ [—a, +a\ is now 

" ' 27ru;K-i(e"« - 1) 27r(i + WK-i) 



F[l,l~iwHr^,2-iwK-^,e-~°"']} (69) 



Here, even in the very low energy regime, the contribution of short distances is not negligible. In fact, it is much 
bigger than the thermal spectrum itself. To see this we can approximate I^^{'wk~^ ^an) as 

I^,{wn \an)^[ + - + 0{{aKY)] - - - 

\TTaK 12-K 7207r J w 2 

an (anY ^ w f {aK^ t'^ \7i 



\ 27r 727r J n \ 72tt J n 

Note that in this case the dominant term is of order l/an, therefore 

\™ ( 2tL-i = — » 1 • 71 

A similar behavior can be found for w ~ Wtypicai- 

We illustrate the difference between both calculations in Fig. 2. With the normal-ordering prescription the short- 
distance contribution is small, in contrast with the ie-prescription. We have chosen a large surface gravity and different 
values of a to better show the effect in the drawings. We clearly observe that, although both prescriptions lead to 
the thermal result when a — > +oo, they do the job in very different ways. 



^ The exponential behavior in frequencies of the ratio I I68I I explains why potential deviations from thermality arise at frequencies much 
lower than to ~ 1/lp. 
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N(w) a=i N(aj) a=10 




Figure 2: Plot comparing the Planckian spectrum N{w,k) = (e^'^""^ — (solid line) with the contributions (dashed 
line) and (dotted line) coming from distances |z| < a according to the normal-ordering prescription and the ie prescription, 
respectively . We have taken k = 0.1 and a = 1, 10, 30 and 100 (in Planck units), respectively. 



B. Fermions 

We shall extend the previous analysis to fermions. The integral involved is 

n+a r (I) 



{wK~^ , an) = — dze~ 

277 



2 



(sinh ^z) 

] u-i + 2«;/^-i)(e-««/2+-- X 

TT 27r(l+ 4^2^-2) 1-'^ ^ '"^ 

13 13 

F\l, iwK~^, iwK~^ ,e^°"^] — F\l, iwK~^, iwK~^,e°"^] x 

L>2 '2 ' ' ^'2 2 



_ga«/2+ia«,^[j^ - + iwK-^, - + iWK'^ , e""])} . (72) 



1 . -1 3 

See Fig. 3 for a graphical representation. Taking into account that aw <C 1 we can expand {wk~^ ,aK) as 

I''{wn-\an) = (M! + 0[(a«)^]) ^ + 0[(a«)^](^)3 + ... (73) 

Note that the term proportional to k/w, appearing in the bosonic case, has disappeared. Therefore, for very low 
frequencies 



< 1 . (74) 



This shows that the contribution of ultrashort distances is negligible, like in the bosonic case. Moreover, for typical 
Hawking frequencies, Wtypical = Th, we have 



I ^{w typical n 

(^g2TTWtypiaalK.^^ _|_ l)~l 



3 • IQ-^ianf «: 1 . (75) 
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This rate is again very small, but the above expressions unravel the fact that the short-distance contribution for 
fermions seems to be smaller than that of bosons. For the latter the contribution of short distances is proportional 
to the first power of na while for fermions it is the third power. 



N(aj) a=l N(aj) a=40 




N(aj) a=1000 N(aj) a=10000 




Figure 3: Plot comparing the Dirac-Fermi distribution (solid line) N{w, k.) — (e^'^""' + with the contribution coming 
from distances |z| < a according to the normal-ordering prescription (dashed line). For completeness we have also plotted the 
result obtained with the ze-prescription (dotted line). We have taken k = 0.1 and a — 1,40, 10^ and 10'' (in Planck units), 
respectively. 



Finally, let us give numerical estimates for relevant astrophysical black holes using the expansion (j73p 

I^{wK-^,aK) _ Q'^wK^(e^^"""'' + 1) 

(g2™«-i + ~ ■ ^ > 

For a black hole of three solar masses, the relative contribution to the total Planckian spectrum is, for w = Wtypicai, 
of order 10^^^'*% and one must go to frequencies of order w/wtypicai ~ 270 to find contributions {wk~^ Jpk) of 
the same order as the total spectrum. For primordial black holes, M ~ 10^^ g, the relative contribution is 10~^^% 
at Wtypicai and we have to reach frequencies of order lA2wtypicai to get a short-distance contribution of order of the 
thermal distribution. In addition to the conclusions stressed in the bosonic case, namely the robustness of Hawking 
thermal radiation for wavelengths of order of the size of the black hole, we have a new result. Fermions seem to be 
less sensitive to ultrashort distance physics than (spinless) bosons. 



V. MODIFYING THE TWO-POINT FUNCTIONS AT SHORT-DISTANCES 



In the previous section, we have investigated the contribution to the Hawking spectrum coming from distances 
z < /p at /+ assuming that the physical laws are not modified at such scales. We found that potential deviations 
from thermality only manifest themselves at high frequencies. One can legitimately wonder, however, why we looked 
at distances at instead of at /~, where the sensitivity of the "in" state to short distances is more apparent. In 
fact, imposing naively a cutoff at /~ has dramatic effects on the radiation due to the enormous redshift caused by 
the horizon (see section |TT]) . We were motivated to impose the cutoff at in order to find agreement with the view 
offered by string theory. The purpose of this section is to shed light on the roles played by distances at and I" by 
using a simple model with a modified two-point function. We shall investigate the potential effects on the radiation 
due to the modified short-distance behavior of the matter field, supposedly coming from unknown physics at the 
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Planck scale. We shall see how our model maintains the robustness of the Hawking thermal spectrum at /+, while 
at the same time being insensitive to sub-Planckian distances at /~ . 

Let us assume that the standard two-point function for the spin zero "in" and "out" vacuum states at /~ and /+, 
respectively, gets modified by new physics at very short distances and becomes 



1 



G° 
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An (ui — ^2)^ + 



(77) 



where a is a parameter of order of the Planck length: a ^ Ip. With this modification the expression for the 
black hole particle production becomes (we omit the transmission coefficients ti{w) and the angular delta functions 
^iii2^mim2 since they are also irrelevant for the discussion of this section) 



HKtl^n) = 4 



1 



dvidv2U^^ {vi)u'. 
1 



An {vi — ^2)^ + 



dv dv 



(78) 



I/- are understood to be the "out" modes and the "out" two-point function, respectively, prop- 
agated back to . Since, according to the standard derivation, the propagation to /~ implies a strong blueshift, 
the modes and G™* might manifest some dependence on the particular details of the modified theory, which 
are unknown to us. Thus, we see no simple way to estimate the form of the modes at I~ . For this reason, it is 
preferable to evaluate the particle production as an integral on /+, as in (|^^ . 



{in\Kt\ 

dvi dv2 
dui du2 



= 4 



I/+ 



/+ 
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An (lii 
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(79) 



I/+ 



is understood to be the "in" two-point function propagated to /+. In this region we can use the standard 



form of the "out" modes u°"*(ui) 



^ ^ since we are considering emission frequencies much lower than the Planck 

/47rto ^ 

frequency wp ~ l/lp- We still have to unravel the evolution of G™ to evaluate the above expression. The modified 
short-distance physics near the horizon could dramatically modify the evolution of the two-point function, so that its 
form at /+ could be rather different from the standard one G™]/- . However, we can make the reasonable assumption 
that the propagation to /+ is affected by new physics in such a way that the short-distance behavior of ^^^^G™ at 
/"'" is identical to that of the two-point function for the "out" state 



lim ^[uAu2)G^\j. 



lim G°"*(mi,M2)|/+ 



(80) 



The above condition can be seen as a natural generalization of the Hadamard condition, i.e., universality of the 
short-distance behavior for all quantum states. The Hadamard condition, which plays a pivotal role in the algebraic 
formulation of QFT in curved spacetime [H , ensures the regularity of expression (j29p to evaluate the Hawking radiation. 
Let us see now how (1801) constraints the evolution of G™ from to /+. Note that G™ can be rewritten as 



Gl: 



G* 



1 + a^Gi' 



(81) 



where G™ is the unmodified two-point function. Since, at late times, G™ evolves according to geometrical optics 
approximation 



J J 1 dill dv2 

G 1/+ = — (wi) — (W2)G I/+ - --^Tl — ^ 1 — v^' 

an an An [v(ui) — v[u2)) 



expression (I5T|) suggests the following evolution for G™ 



dv{u-i) dv{u2) 
du du 



^"l--d^(-)d^;(-)^"l---4^-(.,-.2)^.,«.-.-. 



(82) 



(83) 
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This expression guarantees immediately the Hadamard condition ((80)) . We should stress, however, that the evolution 
of the modified two-point function itself is not equivalent, at least for very small point separations (1*2 — ui)^ ~ a^. 



to the ray tracing (or geometrical optics approximation), which would produce instead (|9T|) (see later) and violate 
the Hadamard condition. For larger separations {u2 — ui)^ ^ the propagation agrees, as it must, with standard 
relativistic field theory and is driven by the large redshift (implying then the usual geometrical optics approximation). 
Plugging this expression in ([79]) we obtain 



-1 



Air'^y/uJ^ Jj+ 



dv{ui) dv{u2) 
du du 



It is worth noting that the modified term 



47rG" 



dvi dv2 



\i+ 



(84) 



(85) 



which can also be regarded as a transformation law under the change v = v{u), guarantees the absence of particle 
production under the same group of symmetry transformations (Mobius rescalings) as those of the theory with 
a = 0.^° Assuming that the geometry remains classical (the black hole scale k is well above the Planck scale a), we 
can use in ([M]) the expression v{u) = vh~ K^^e^''", which represents the relation between the "in" and "out" inertial 
coordinates. Performing then the integration in U2 + ui, we are left with {z = U2 — ui) 



{in\N° 



1 



\in) 



2ttw 



dze 



1 



(sinh |z)2 



(f)2a2 



+ a 



Finally, performing the integration in the complex plane, the particle production rate becomes 



{^n\Ktn\^n) = 
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where 



= arctan 



any^l — a^K^ /A 
(l-a2^V2) 



(86) 
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(88) 



The thermal Planckian spectrum is smoothly recovered in the limit a — > 0. Moreover, for a ^ Ip, the deviation to 
the thermal spectrum is negligible for small values of wk~^. This deviation can be expanded as 



2'KWK 



1) 



l)-i 



(89) 



For astrophysical black holes, Ka <C 1, the second factor is negligible for frequencies up to ^ '^G'^w typical, in complete 
agreement with the results obtained in section Hvl (compare, for instance, with (f68|) ). 

The above discussion shows that, despite the apparent sensitivity of Hawking radiation to high energy physics (see 
section 111)) . a Planck-scale modification of the two-point function does not necessarily imply a substantial change of 
the Planckian spectrum. This is so, at least, if the modified two-point function obeys a modified Hadamard-type 
condition. The simplest realization of this condition turns out to be equivalent to the preservation of the powerful 
conformal (Mobius) symmetry existing in the unmodified theory. This seems an unavoidable requirement if the 
corrections to the Planckian spectrum are to be in agreement with the results of string theory in the low-frequency 
limit u; — > 0. The effect of the generalized Hadamard condition is to constrain the short-distance behavior of the 
propagated "in" two-point function, G™|/+, in such a way that it remains close to — l/(47ra2) through its evolution 



Even more, it is what exactly guarantees the invariance of the production rate, up to a shift on the emission frequency, under a radial 
boost with rapidity u — ► u = e^u, v ^ v = e~^v. 
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to /"'", despite the large blueshift. In fact, G™|/+ is an observer-independent quantity in the limit xi X2, i.e., it 
tends to — l/(47ra^) for any function v = v{u). Note in passing that this condition is somewhat related to approaches 
to quantum gravity aimed at deforming Lorentz symmetry while keeping the principle of relativity [l^ . 

To conclude, we note that if the deformed two-point function at 

An [vi — V2) + 

is naively propagated (i.e., by ray tracing) to I"*" as 

-J dvi dv2 

Qi^^j^ I c^"! du2 (91) 

" Att {v{ui) — v{u2))^ + ' 

where (5™!/+ = ^(ui)^(it2)G™|/+, the particle production rate is now time-dependent and the thermal spectrum 
is lost for any nonvanishing a. 



VI. CONCLUSIONS AND FINAL COMMENTS 

It is highly non-trivial [H to truncate Hawking's derivation of black hole radiance to account for unknown physics 
at the Planck scale. A simple estimate of the contribution of virtual high frequencies apparently shows that they 
are essential to produce the thermal outcome. One can then change strategy and try to evaluate the contribution of 
Planckian physics in position space, which requires a rederivation of the Hawking calculation in terms of two-point 
functions, as we have explicitly shown in section [TTll When these two-point functions are treated in the distributional 
sense, with the usual ie prescription, one reproduces exactly the thermal result. However, one can equivalently handle 
the divergence of the two-point function by trivially taking normal ordering. The consistency of this procedure 
is guaranteed by the Hadamard condition: the short-distance behavior is universal for all physical states. The 
advantage of this second option is that it offers a natural way to evaluate the contribution of short distances at 
to Hawking radiation. 

We have found that the contribution of short-distances at low frequencies w <^ k is negligible. Our analysis allows 
us to go further and investigate the short-distance contribution for frequencies of order the Hawking temperature 
Th and beyond. We find that the contribution of ultrashort distances is also negligible for frequencies of order Th- 
In fact, for a black hole of three solar masses we need to look at high frequencies, w/wtypicai ~ 96 (for bosons) or 
w / Wtypicai ~ 270 (for fermions), to find that the contribution of Planck distances is of order of the total spectrum 
itself. This suggests that Hawking thermal radiation is very robust, as it has been confirmed in completely different 
analyses based on black hole analogues; in string theory (for large wavelength) in near-extremal charged black holes; 
and also in some models of canonical quantum gravity [23| . 

One can legitimately ask why, in section IV, we evaluate distances at instead of just at where the sensitivity 
of the "in" state to high energy scales is more apparent, as we showed in section |TT1 Our heuristic motivation is based 
on the view offered by string theory, where the Hawking radiation is obtained as the result of collisions between open 
string excitations. In that approach, the standard large blueshift of low-energy gravity theory does not seem to play 
the pivotal role that it does in the pure semiclassical treatment. The fact that we consider the fundamental Planck 
scale at does not immediately guarantee that the Hawking radiation is kept unaltered from Planck-scale physics. 
As we show in section IV, with the standard ze-prescription the short-distance contribution to Hawking radiation is 
not negligible. In contrast, with the normal-ordering prescription the bulk of the Hawking effect is maintained at low 
frequencies, in agreement with the results of string theory. 

In addition to the above arguments we have approached the problem in section V in a different way. We have 
considered an explicit modification of the two-point function at the Planck scale. Motivated by the crucial role played 
by the Hadamard condition in the ordinary relativistic theory, we have assumed that the short-distance behavior of 
the modified theory should also satisfy a sort of generalized Hadamard condition (universal short-distance behavior). 
The simplest realization of this idea turns out to be equivalent to the preservation of the powerful conformal (Mobius) 
symmetry existing in the unmodified theory. Armed with this condition, the contribution to the particle production 
rate of the "in" and "out" two-point functions in (|84p is similar when they are compared in the same ultrashort range of 
distances, despite the large blueshift horizon effect. As a result, the two contributions compensate each other and lead 
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to an emission spectrum very insensitive to trans-Planckian physics. The generahzed Hadamard condition, therefore, 
seems to be necessary to maintain the bulk of the Hawking effect. Moreover, it is in this context that the apparent ten- 
sion between /+ and to measure separations is elhminated since in both we find the same finite short distance hmit. 

A last comment is now in order. In the string theory analysis one has, at least, two relevant parameters: the surface 
gravity k and the radius Vg of the supersymmetric, charged black hole. The surface gravity is assumed to be small, in 
comparison with the inverse of the size of the black hole, i.e., k <C l/?'g. The emission frequency can reach n, but can 
never reach l/vg (or become larger) to guarantee the validity of the string theory calculation. Obviously the analysis 
of string theory excludes astrophysical black holes of the Schwarzschild type (for which k ^ ^/fg)- Our results, 
however, suggest that one could also expect string theory to predict in this case, in some subtle way, agreement with 
Hawking's results for frequencies around l/vg and, at least, a few orders beyond. This is so because we do not observe 
any significant contribution to the thermal spectrum coming from the short-distance region, where new physics could 
arise, up to such high frequencies. This fact offers a very non-trivial challenge for any quantum theory of gravity 
having computational rules very different from those of semiclassical gravity (as in string theory or background- 
independent approaches), since when w ~ l/j'g the grey-body factors Ti['w) cannot be computed analytically. They 
are only known numerically |21|] , as a result of solving field wave equations in the black hole background. String theory 
manages to account for the greybody factors in the low-energy regime, where they admit an analytic expression. In 
fact, for all spherically symmetric black holes, the low-energy absorption cross section is proportional to the area 
of the horizon 22, 23]. But for typical Hawking frequencies the grey-body factors remain elusive for any analytic 
treatment. Reobtaining them from such a different computation would be extremely impressive. 



Appendix A 

We will complete here the steps missing in the derivation that led to the emission rate (P5|) . Using the wave packets 
we can express it as 

(^'^l^iinijsnsl*'^) = / dw' Pj,n,,w'P%n2,w' = (92) 
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Using (|20)) we get 



p-i{wi-W2)vH _ _ 

r(l + iKr^w{)V{\ - in-^W2)Sa\Kr^{wx - W2)\ . (93) 

This integral can be estimated explicitly when the width e of the frequency interval [je, {j + l)e] is assumed, as usual, 
small. In this case, the integral is essentially as follows 



where 

e/2 i.e/2 



4i«2('^) = -/ dxi d^^^^[^-v„]x^-i[^-VH]x2-^K (^i+^2)/2^^[^-l(^^ „^2)] (95) 

-e/2 J-c/2 



and xi^2 = wi^2 — {j + l/2)e. The factor Sj-^j^ in ([M|) is due to the role of Sa, which selects frequencies on a very 
narrow band of order \wi — W2\ ^ na. For this reason, it is also convenient to introduce a new variable y — xi — X2 
and rewrite /„j„2((t) as follows 

Imn2{<y)~- dxie / dye'[^-^'«l^<5,[K~iy] . (96) 

C J-e/2 Jxi-e/2 
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In writing this we have neglected the term e"'^ {xi+x2)/2 t^jY^-^^^ jg almost constant (unity) over the integral. We can 
now estimate the integral over y having in mind that Jo- is very well approximated by a square step of width ttkct and 
height l/(7r cr) centered at y = 0. This means that the main contribution comes from the interval [— ^^]. This 
fact makes the outcome of the integral independent of .ti, which also allows us to perform the integral in xi. Putting 
all together we find 

Plugging this result back into ([94]) we find ([25| . 

We will now briefly consider the effect of introducing the cutoff in frequencies in a different way. The cutoff was 
introduced in ((20l) in the form 

^log[A/K] 
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We will now consider the change 
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where A must be of order ^ log[A/K]. This modification leads to a redefinition of 
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5^[k-\wi - W2)] = ^ — ^ (100) 

which in the limit 2a ^ Q also becomes Dirac's delta function. One can then proceed as above and define the 
corresponding function I„-^n2{^)y which this time can be evaluated extending up to infinity the limits of integration 
over the variable y — xi ~ X2- This leads to 

/,u„,(a) = At(5„,„,e-[(^-''«)^*]' (101) 

The corresponding emission rate is now 

(mliV""*'" \in)=S- S _IM!!!lll!_e-[(^-"«)'^'*]' (102) 



This expression is always positive definite and exhibits the same decay rate as p5p if we identify 2a with a, which in 
fact is the right choice for the definition of (jlOOp . 



Appendix B 

We will proceed now to solve the massless Dirac equation in a curved background with spherical symmetry. The 
equation to solve is^^ 

7^V^V = (103) 

where 7^ = j^VJ^ix) satisfy {7'', 7"} = 2^^"^, {7^,7''} = 2ry"'' and VJ^Vj^^ri"'' = g^"" represent the vierheins^'^ . Note 
that V^iV = (c^M ~ r^)'!/' where T^ = —j'^^^'^V^^^V^c represents the spin connection. We will take the curved space 



For earlier references see [24|. and for a more advanced treatment (no needed for the purposed of this paper) see [25ll . 
Due to our convention for the metric signature the 7" matrices should verify the conditions (7")^ = — /, (7')^ = I- However, to agree 
with the standard notation for Dirac matrices in this appendix we have flipped the metric signature to (+, — , — , — ). This is, however, 
irrelevant for the computations carried out in the body of this paper. 
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line element ds"^ = e^Pdx'^dx — r^dft^, with dx~^dx = dt^ — dr*"^, and rjab = diag{l,— 1,-1,-1). Introducing the 
ansatz ip — $ and making the simplest choice for vierbeins (i.e., to be parallel to the unit vectors in t,r*,9,(j) 
directions), the Dirac equation (|103p becomes 
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where the index i runs over the non-angular variables. Since for x^ = t zt r* we have j'^V^di — e ^[y^dt + 7^9^*], 
(|104p can be written in the more familiar form 
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The angular part of this equation can be reexpressed as eP^^K/r 
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where the operator K 
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commutes with the Dirac equation as well as and J3 and, therefore, its eigenvalues can be used to characterize the 
angular part Xm^Kj of the modes: KxnijKj = (-'tjOXm^Kj, with k| = (j + i)2^ Moreover the eigenfunctions Xm^K^. 
admit the following decomposition XrnjKj — c^XmjKj + ^ Xm^K^i with 



(105) 



(106) 



(107) 



'7(^)"' 








Therefore, in a stationary spacetime, p = p{r), a general solution can then be expressed as 



(108) 
(109) 



-iwt 



G„.,(r)^(r)™/ 



(110) 



where we have used that a^r]{f)^. — r]{r)^^^ . and the functions F^Kj (f^*) and Gw^j (f^*) satisfy the following equations 
(see also [23| ) 



(111) 

(112) 



Adding the time-dependent part, the above equations lead to plane-wave solutions ^ e ) = e ™^ for all kj 

as r — > 00. 

We note that the form of the eigenfunctions XmjKj can be worked out immediately if the vierbeins are chosen to 
be parallel to unit vectors in the standard t,x,y,z directions''^. The bispinors rjmjKj can be constructed, as it is 



With this orientation for the vierbeins K can be written as if = 7" (7 -|- 25 ■ L) , which is the standard form of this operator in Minkowski 
space. 
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well-known, using the Clebsch- Gordon rules for addition of angular momentum in terms of spherical harmonics and 
two-component spinors, and the result is 



and 



j+rrij -ymj-1/2. 



2j 



i-1/2 



(113) 



2j+2 -^3 + 1/2 \^ 
j+l+nij ymj + l/2, 
2j+2 -^3-1-1/2 



(114) 



With the modes given in (jllOp conveniently normalized, the quantized Dirac field can be expanded in modes as 

^^h. L,.,„,„„.,„,(.) + .L,„,...,.»l (115) 



where Uwn m-ix) and VwKm {x) represent positive and negative energy solutions respectively. On the other 
hand, since we are dealing with massless spinors, it is necessary, on physical grounds, to use states with well 
defined helicity. In particular, left-handed spinors can be obtained from (jllOp by projecting with = i(/ — 7^), 
where 7^ = i7*'7-'^7^7^. We will therefore be working with the (normalized) modes ipwjmj = '^{'^w\Kj\mj —'4'w~\Kj\mj)- 

We will now carry out the calculations that lead to ((52|) (adapted now for chiral spinors). First thing to note is 
that the propagated backwards mode (jH]) contains a term of the form yjdu{v)/dv. A simple way to realize why this 
term arises is that it is necessary to ensure the invariance of the scalar product under time evolution. Putting aside 
backscattering effects, the Dirac scalar product for out modes can be written, equivalently, as 

dr2dur2u™*7+u°"* = / dnr^dvu°''*j-u'"'* . (116) 
/+ Ji- 

The above equality requires, up to relative signs in the spinor components, that 



u'"'\v)\j+ = ^/du{v)/dve{vH - 'y)u°"*(M)|/- . (117) 

Note that the factor e^P^^ in ()110|) also signals this behavior. Since the spinor ^{x) does behave as a scalar under 
general changes of coordinates, it follows that the functions F and G must somehow compensate the change in e"''/^ 
under conformal transformations. 

Let us now focus on the integration over the angular variables prior to (j52p . This integration can be readily 
performed if we put the result of ([501) i^i^o (|T7)) . We then find 

X dv.rldn, i^vir^^(:r,)^jl-2^ur^'^''ix,)^ (118) 
where the indices Ji,i2 and k denote {w,j,mj). Using the modes of eas. (|49p and (|5ip it is immediate to verify that 



-out,Lf J7°-7^] in,L( . ^32(^2) du{v) 

dn^U^^ ix2) Vk' {X2) = ^^V^^^"''""^ 

X e''"'"^^''>+''"^'S,n,.^,nJnj, (119) 

where we have used that J dfi r]"^/\r)r]^^^,' (f) ~ Smjm.j,SKjK^, ■ An analogous calculation applies to the second factor 
in (jllSp . Plugging these results back into (|118p we obtain 

47r2 



{in\Ni^^^\in) = Smj^m,^Jjij2— j^r / dvidv2 (120) 



V dv dv Jq 
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There remains to perform the integration in w, which yields 







e-*o [vi -V2- ie) 



and leads to the sought-after result. 
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